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We propose a possibility of ultra-relativistic electromagnetic counterparts to gravitational waves 
from binary neutron star mergers at any viewing angle. Mechanisms include the merger-shock 
propagation accelerating a smaller fraction of the neutron star surface to larger Lorentz factor, with 
the outer parts carrying less energy, ~ 10 47 r _1 erg. This mechanism is difficult to resolve by current 
3D numerical simulations, and depends on the neutron star equation of state. The outflows emit 
synchrotron flares for seconds to days by shocking the surrounding medium. Ultra-relativistic flares 
shine at early times and high energies, potentially detectable by current X-ray to radio instruments, 
such as Swift XRT and Pan-STARRS, and even in low ambient density ~ 10~ 2 cm" 3 by EVLA. The 
flares probe the merger position and time, and the merger types as black hole-neutron star outflows 
would be non/mildly-relativistic. 

PACS numbers: 04.30.Tv,97.60. Jd 



Introduction. Binary neutron star (BNS) mergers are 
main sources of gravitational waves (GWs) for ground- 
based laser-interferometric detectors, such as advanced 
LIGO, advanced Virgo, and KAGRA (formerly LCGT), 
in the coming five years PHI]- GWs test gravity, open a 
new window for astronomy, and probe the supranuclear- 
density matter in neutron stars (NSs). We expect a few 
tens events per year within a few hundred Mpc 

A simultaneous detection of electromagnetic (EM) sig- 
nals is indispensable for declaring a confident discovery 
of GWs Since "hearing" a sound of GW entails 

a bad localization about degrees at best, "seeing " the 
EM counterparts will also not only increase the GW sen- 
sitivity but also expand multi-messenger astronomy by 
finding the host galaxy, its redshift, etc. 

Short gamma-ray bursts (SGRBs) are plausible coun- 
terparts to BNS mergers [7] . GWs will verify the merger 
hypothesis for SGRBs. However, SGRBs are usually 
strongly beamed into a small angle Q. Most SGRBs 
are off-axis to us without detectable emission, although 
GW observation is biased toward the binary's rotational 
axis (i.e., probably the jet axis). The "orphan" afterglow 
is also dim, which is produced by the off-axis jet deceler- 
ated down to a mildly relativistic speed at late time. 

Two promising models have been proposed for nearly- 
isotropic EM counterparts 0, Both of them are 
based on non/mildly-relativistic outflows with ~ 0.2-0. 3c 
(roughly an escape velocity of the NS) from a compact bi- 
nary merger. The outflows can be produced by neutrino- 
driven wind , magnetically-driven wind [ll], [l2[ , tidal 
ejection and shock- wave ejection [3] (see also be- 
low). Recent fully general relativistic simulations show 
the ejecta mass of > 10~ 3 M Q for a wide range of pa- 
rameters even without neutrino nor magnetic effects [15] . 
One emission model is the macronova (or kilonova) which 
shines on ~ day after the merger in UV-optical bands 



via radioactive decay of r-process elements 

[Mi. The 

other model is radio synchrotron emission from the col- 
lisions between the ejecta and the ambient medium, like 
GRB afterglows, after ~ year from the merger [1, 

In this Letter, we suggest a possibility of nearly omni- 
directional ultra-relativistic counterparts to BNS merg- 
ers, even if the GRB jet is off-axis to us (for the first 
time to our knowledge). We consider the shock waves 
produced by the BNS collision [l!| |2(|. Shock waves 
propagate in the NS non-relativistically at first, but ac- 
celerate to a relativistic velocity down a steep density 
gradient at the NS surface. After the shock breakout 
from the surface, the shocked matter expands into nearly 
a vacuum, converting internal energy into kinetic energy. 
The resulting Lorentz factor T is larger for outer and less 
energetic parts of ejecta. Such a trans-relativistic accel- 
eration has been discussed in the context of supernovae 

We estimate the relativistic ejecta mass to be ~ 
10~ 7 r~ 2 Mq for r » 1, and calculate the synchrotron 
radiation from a relativistic blast wave decelerated by 
the ambient medium and energized progressively by the 
inner catching-up outflows with a range of T. The ultra- 
relativistic nature makes the flare bright at early times 
(from second) and high energies (to X-ray) in contrast 
to the non/mildly-relativistic cases. We find flares are 
detectable by current X-ray, optical, and radio observa- 
tories, such as Swift XRT, Pan-STARRS, and EVLA for 
our fiducial case. 

The counterparts signal the merger time more precisely 
than non/mildly-relativistic models. The counterparts 
could also enable us to distinguish the merger types be- 
cause black hole-NS mergers are not likely to accompany 
strong shocks for ultra-relativistic outflows. 

A point to be noted is that the 3D numerical sim- 
ulations of BNS mergers have not sufficiently resolved 
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the NS surface. Although the results for non/mildly- 
relativistic ejecta are solid and the existence of shock 
waves is implied by the heatup of the colliding region 
[l9l l20j |. currently it is not feasible to follow a small 
surface mass to a ultra-relativistic speed because of the 
numerical viscosity, artificial atmosphere, and limited 
computational resources, even for Newtonian gravity. 
Therefore, it is worthwhile to highlight the impacts of 
ultra-relativistic outflows for motivating the future well- 
resolved calculations. 

Acceleration. We first consider the mass ejection just 
after the NS collision. The NSs come close as a result of 
the GW radiation reaction, and collide with each other. 
The colliding part is shock-heated up to a temperature 
fcfiT ~ 50 MeV. Because of the oblique collision, the 
shocked region has a pancake-like shape with its thickness 
R s h ~ O(l) km [H, Ho|, as shown in Fig. [TJ 



The hot matter tries to expand to a cold, low-pressure 
region, i.e., to the NS surface not to the NS core (see the 
top panel of Fig. [1]). This is like a tomato that splashes 
sideways when dropped on the floor. The shock-wave 
velocity is comparable to the sound velocity of the core 
material, Ui n ; ~ 0.25c [25j. At this stage, the shocked 
matter cannot escape from the merged BNS, because the 
expanding velocity is less than the escape velocity, 
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where M» and i?* are the typical mass and radius of the 
merged BNS. 

The shock is accelerated down a steep density gra- 
dient of the NS surface. The shock velocity increases 
as oc p~ Qnl with cv nr ~ 0.187 in the non-relativistic 
regime, where the index depends weakly on the equa- 
tion of state (EOS) [2l|. Once the shock is accelerated 
beyond ~ 0.5i> CS c7 the shocked material can escape from 
the BNS by converting thermal energy into kinetic energy 
to obtain « v csc later 2l|, 23 1. 



The shock is initially weak, because the shock velocity 
is comparable to the sound velocity, but becomes strong 
as it propagates in the NS crust where the sound veloc- 
ity decreases outward [25| . The density profile of the NS 
crust is approximately given by p oc x™, where x is the 
depth from the NS surface (p = 0) and n is the poly- 
tropic index of the NS crust matter. We adopt n = 3 for 
a fiducial value, because this is for the relativistic degen- 
erate electron gas and is consistent with more detailed 
nuclear-theory-based EOSs of the NS crust [26j |. 

The shock velocity increases by a factor of 
0.5i>esc/wini ~ 1-5 when the density drops by l^" 1 /" 1 " 7 « 
0.1. The mass outside the corresponding density can es- 
cape as ejecta. For the density profile p oc x n , the ejecta 
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FIG. 1. Top: A schematic picture of the shock generation, 
propagation, and mass ejection at onset of the BNS merger. 
Two shaded (blue) ellipses are the BNS with a low tempera- 
ture, and the low-density crust is depicted with fainter shade 
(light blue). A dark (red) region at the contact surface is the 
shock-heated region. Black arrows denote the BNS motion 
just before the merger. Shock waves are generated from the 
contact of the BNS. The shocks become strong near the NS 
surface and eject a part of the NS crust ultra-relativistically. 
Bottom: A snapshot of merging BNS with 1.5 Mq taken from 
a simulation in The temperature on the equatorial plane 
is shown, and the contact surface is heated up to ~ 50 MeV. 



mass is estimated as 

R* V Wini 



-7.1 



-"xlO-i^Jgj)"^)- , (2) 
where ~ 5.35, and a factor i? s h/i?* ~ 1 km/15 km 
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FIG. 2. The kinetic energy distribution of ejecta with a ve- 
locity larger than fSF for various polytropic indices n and 
the ejected mass by the shock breakout M a h- We assume 
Af* = 2.8M© and R, = 15 km. 



is the geometrical fraction of the crust mass M crus t ~ 
O.OIM0 that is swept by the shock (see Fig. [T]). 

The outer, less massive ejecta is faster up to a ultra- 
relativistic velocity. For example, the shock attains Y s ~ 

10 as the density drops by a factor of ~ Y^ +i ^ 3 ^ iq 4 
within a thin layer ~ -R C rust/ri 2+4/v ^ )/n ~ i? crust /30 
[jjj ]. where i? C rust ~ 1 km is the crust thickness. Then, 
after the breakout from the NS surface, the shocked mat- 
ter accelerates to a Lorentz factor T ~ r 4+v/ ^ > 500 by 
converting heat into motion with the aid of the pressure 
of the inner ejecta. In order to resolve the thin layer 
< -Rcrust/30, a grid size of < 10 m is required for the nu- 
merical simulations. Such a high resolution is unrealistic 
for current and near-future BNS simulations in 3D. 

To make a detailed estimate, we apply the trans- 
relativistic acceleration of a supernova exploding the stel- 
lar envelope [24j. The kinetic energy of ejecta with a 
velocity above (5Y is given by Eq. (56) of [24|, 



E{> PT) = pSS" F{pT)M sh c 2 



(3) 



Here, p csc is the shock intensity required for the shocked 
material to escape, and is evaluated by /3 CSC w 1.73p C sc 
27j according to_24|. The distribution F(/3T) is given 



by Eq. (38) of [24[. The gravitational binding energy of 
the merged BNS does not contribute significantly to the 
ultra- relativistic component. 

Figure [D shows the kinetic energy distribution of the 
ejecta E(> f3Y). We can see the ultra-relativistic energy 
E(> fir » T = 10) > 10 46 erg for our fiducial case. For 
/3T > 1, Eq. © yields 



E(>T) « 2.6 x 10 47 
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where we assume M* = 2.8M Q , i?* = 15 km and n = 3. 
The corresponding mass is ~ 

10 -7 r -2.i MQ for T ^> 1. 

The high-r component carries small but still appreciable 
energy for the emission as E{> Y) cx r -0 58-1 - 58 /™ ~ 
T~ . The energy distribution becomes harder for larger 
n, providing a possible way to infer the EOS of the NS 
crust. 

The energy distribution is sensitive to the value of 
Af s h, and therefore Di n i and the polytropic index, n. 
The density profile could also be affected by the neu- 
trino/magnetic wind and the tidal torque. If the density 
at the surface does not go to zero rapidly as in the stellar- 
wind density distribution, the shock acceleration will not 
work efficiently. The break down of plane-parallel ap- 
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proximation in 24| could modify the Y distribution. The 
ultra-relativistic ejecta may be decelerated before emis- 
sion by the surrounding material, such as tidally-ejected 
material before the merger. The accurate estimation of 
the ejecta mass M s h and the Y distribution is left as the 
future study. 

The ejecta may be also accelerated when the density 
waves propagate across the entire core to the opposite 
surface. The amount of shock-breakout ejection would 
be larger by an order of magnitude than Eq. ([2]) since the 
geometrical fraction is large. As a result, the emission is 
brighter and longer. Whether density waves can propa- 
gate across the core or not depends on the NS structure, 
and thus on the EOS of the NS core. 

Radiation. Next, we calculate the spectra and light 
curves of EM signals applying the synchrotron shock 
model of the GRB afterglow 28, 2(|. The outflow gener- 
ates a forward shock sweeping the surrounding medium 
with a constant number density tih- A fraction es ~ 0.01 
of the internal energy released by the shock amplifies the 
magnetic field, while a fraction e e ~ 0.1 accelerates elec- 
trons with a Lorentz factor distribution dN e /d^ e oc 7~ p 
for 7 e > j m and p ss 2.2, which emit the synchrotron 
radiation. 

We assume that the ejecta outspreads in a nearly 
spherical manner after the breakout. If the ejecta is con- 
centrated within an angle 6 from the annular shocked re- 
gion (see Fig.Q]), an observable region decreases by « 0/w 
but the isotropic energy increases by ~ tt/9 for a single 
BNS merger, and finally detection rates will increase by 

We take 100 Mpc for a fiducial distance to a BNS, given 
the plausible merger rate ~ 40yr~ 1 within the angular- 
averaged horizon ~ 200 Mpc [4|. We neglect the cosmo- 
logical effects, since the maximum horizon is 445 Mpc for 
advanced GW detectors. Our fiducial model can be de- 
tected up to > 200 Mpc for optimal observational bands 
and 7Jh as shown later. 

The outflow carries larger energy in lower-T (inner) 
part in Eqs. ((3]) and (j4|. The slow flow rear-ends and 
refreshes the external shock which is decelerated by the 
ambient density [3(J 31 1. We illustrate formulae for a 
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single power-law distribution of the kinetic energy, 

E(>T) = EY 1 "*, (5) 

where E = 2.6 x 10 47 erg and s = 2.1 for our fidu- 
cial case with r > 1, and E and s rise as V drops in 
Eq. ([JJ. We also adopt an ultra-relativistic approxima- 
tion, r > 1. Once the fastest flow begins to decelerate, 
the catching-up condition for a slower flow is given by 
ET 1 - 3 = 167rr 2 -R 3 n H mpC 2 /17 [El, HI]. Since the shock 
radius is connected with the observer time as R « 4Y 2 ct 
[28j . we obtain the hydrodynamic evolutions 

(ct/ls)- 3 ^, 
4Y(t) 2 ct, 



r(i) 

R(t) 



(6) 
(7) 



,2)1/3 

Note T oc t~"' si and R oc i 034 for our fiducial s = 2.f, 
compared to T oc i -3 / 8 and R oc t 1 / 4 for a single- velocity 
outflow. 

Given the hydrodynamics model above, we can calcu- 
late the evolution of radiation 



where i s = (HE /lQ2A-KnY{m p c 2 ) lri is the Sedov length. 
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29( | . The synchrotron 
flux has a broken power-law spectrum, F v oc zW 3 for 



v < v m and F,. 



/-0-l)/2 



for f > z/™. The 



cooling frequency is high (>MeV) for typical parame- 
ters. The self- absorption frequency is only relevant for 
the radio band with high hh ^ 1 cm -3 . The characteris- 
tic synchrotron frequency v m oc t _12 /( ;5 + 7 ) decreases with 
time. As v m crosses the observed frequency v, the flux 
reaches the maximum value with the peak time and flux 
given by 



,1.52 ,0.38 p 1 / 3 ^ 0.05, -0.76 



£peak = 6.2 mS ^ e -i^B, -2-^47 n H "18 

peak — U.^l ^iJy e e _ x e B _ 2 &47n H u 2 v 1& 



(8) 
(9) 



for fiducial parameters, where Q x = Q/10 X in units 
of erg for E, Mpc for D, and Hz for v. The flux 
grows as F v oc i( 3s+1 )/( s + 7 ) w t - 80 and decays as 



F„ 



oc 
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-0.43 



across the peak. These 



light curve behavior could constrain the T distribution, 
i.e., s, and hence the crust EOS, n ps 1.58/(s — 1.58), 
in principle, although the Y distribution may be modi- 
fied by the neglected effects, such as the breakdown of 
the plane-parallel approximation. The degeneracy be- 
tween n and p is solved by the high-energy spectral index, 



oc v 



-(p-i)/2 



Figure [3] shows the light curves in X-ray, optical, and 
radio bands using the (non power-law) velocity distribu- 
tion in Eqs. ((3|) and (j4j with various ISM densities 71h- 
In contrast to the non/mildly-relativistic cases, i.e., op- 
tical macronovae and radio flares, the ultra-relativistic 
signals appear in the early phase down to seconds and in 
the high frequency up to X-ray. The X-ray and optical 
peaks correspond to Y m 400 and 70, respectively. The 
flux becomes smaller for lower ISM density. We also plot 
the analytic light curves of the single power-law distribu- 
tion in Eq. ([3]) with nn = 1 cm -3 , which underestimates 
the late (low-r) flux. 
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FIG. 3. Light curves of the ultra-relativistic EM counter- 
part to the BNS merger with E = 2.6 x 10 47 erg at 100 
Mpc distance in the X-ray (0.2-10 keV integrated flux in 
erg cm -2 s -1 ), optical (in 629 nm, r-band magnitude), and 
radio (1.4GHz in fiJy) bands for various values of nu to 
r ~ 1. Long-dashed (red) curves are the analytic approx- 
imations with the single power-law distribution in Eq. © 
with nn = Iran" 3 . Short-dashed (black) lines are sensi- 
tivity curves of Swift XRT, Pan-STARRS, LSST, ASKAP, 
and EVLA. Dot-dotted (black) curves in the middle and bot- 
tom panels show the non-relativistic EM counterparts of the 
macronova in bolometric magnitude and the radio flare [f| , re- 
spectively, with the ejecta mass of 1O" 3 M , velocity /? cj = 0.3 
and «h = 1 cm -3 . 



Discussion. The ultra-relativistic counterpart pro- 
posed here is bright from the early epoch soon after the 
BNS merger, and decays rapidly. While this feature is 
advantageous to confirm a tight association with GWs, 
the observation will require effective strategies. One such 
strategy is full-time EM monitoring of nearby (up to 
~ 100 Mpc) galaxies, where the EM signals trigger the 
GW analyses like the SGRB case. The method is expen- 
sive but find many other transients including supernovae 
as a by-product. 

The other strategy is prompt follow-up by EM instru- 
ments in response to rapid alerts from GW detector net- 
works. The localization requires at least three and hope- 
fully more than four GW detectors. Because a localiza- 
tion error will be ~ 1 degree 2 at best for a BNS merger 
33| , covering this large area is crucial for the EM follow- 
up. SwiftXRT has (0.4 degree) 2 FOV, and can detect the 
decaying phase of the X-ray signal by tiling the FOV. De- 
tecting the X-ray peak may be possible if EM precursors 
are observed in advance 34| . The optical flare can be ob- 
served around its peak by all-sky surveys, such as Pan- 
Starrs with 7 degree 2 FOV and LSST with 9.6 degree 2 



FOV, if nn > 10 cm . The radio flare can be also 
detected around the peak by EVLA even in low ambi- 
ent density nn ~ 10~ 2 cm~ 3 . The low ambient de nsity 



is suggested by the radio observations of SGRBs [35 
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Since EVLA has relatively small FOV of 0.25 degree 2 , 
large FOV instruments such as ASKAP may be more 
realistic choices. 

We also speculate that GeV-TeV 7-rays could be gen- 
erated via inverse Compton scatterings or hadronic pro- 
cesses, such as p-7 collisions. If energy E 1 is converted 
to 7-rays with typical energy e 7 , an expected number of 



photons N 1 
be 



for a detector with area A on the earth will 
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This suggests that the km future instrument such as 
CTA could also detect the EM signal in 7-rays. The TeV 
7-rays are not attenuated by the infrared background at 
~ 100 Mpc. The other multi-messenger signals such as 
neutrinos and cosmic-rays may be also worth to examine 
in the future. 

Ultra-relativistic outflows could also be produced by 
other mechanisms such as the Poynting wind from the 
NS surface with a small baryon load like in the magne- 
tar models for GRBs [36j |. In this case the EM signals 
could arise from the magnetic reconnection without the 
radioactivity or the ambient medium. 
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